Antiferromagnetic spin ladders effectively coupled by one-dimensional electron liquids 
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We study a model of the stripe state in strongly correlated systems consisting of an array of 
antiferromagnetic spin ladders, each with ni eg legs, coupled to each other through the spin-exchange 
interaction to charged stripes in between each pair of ladders. The charged stripes are assumed 
to be Luttinger liquids in a spin-gap regime (Luther- Emery) . An effective interaction for a pair of 
neighboring ladders is calculated by integrating out the gapped spin degree of freedom in the charged 
stripe. The low energy effective theory of each ladder is the usual nonlinear <r-model with additional 
cross couplings of neighboring ladders. These interactions are found to favor either in-phase or 
anti-phase short range spin orderings depending on whether the charge stripe is site-centered or 
bond-centered as well as on its filling factor and other physical parameters of the charged stripe. 



Experiments in high T c cuprates and other doped 
antiferromagnets have established that electronic stripe 
ordering exists in these materials [0. In this phase, 
the doped holes are confined into quasi-one-dimensional 
structures separating locally antiferromagnetic (AF) re- 
gions. In addition, there is by now extensive evidence 
suggesting that dynamical or static stripe order persist 
in both the normal and superconducting states of the un- 
derdoped cuprates such as La2- x Sr x Cu0 4 . @ These ex- 
perimental findings have led support to the picture that 
stripe order may be intimately related to the mechanism 
of high T c superconductivity. For this, it is important to 
understand the mechanisms of stripe formation. 

The interplay between the charge stripe and the in- 
tervening Mott-insulating regions has received some at- 
tention. The spin-gap proximity effect was found to be 
induced by pair hopping between the stripe and its en- 
vironment in the theory of Emery, Kivelson and Zachar. 
H A Landau theory was also proposed by Zachar, Kivel- 
son, and Emery Q] to describe the intermediate rela- 
tion of the charge stripe order and the incommensurate 
spin modulation. They argued that the stripe phase is 
charge driven through continuous or first-order phase 
transitions: charge order sets in first and then anti- 
phase spin domains set later. Their interpretation qual- 
itatively agrees with the experimental observations of 
Lai.6-a;Ndo.4Sr x Cu04 Q and La2NiC>4+<5 [[§■ Granath 
and Johannesson [Q studied a Hubbard chain (represent- 
ing the stripe) coupled to an AF spin ladder and found 
that the magnetic correlations on the ladder generate a 
spin gap on the stripe. Very recently, Zachar Q pro- 
posed a domain wall (stripe) model of electron dynamics 
in the AF environment to investigate the energetics of 
anti-phase and in-phase spin modulation. 

In this paper, we attempt to find the effective inter- 
action between two neighboring spin regions mediated 
by the charge stripe. Such a calculation is important 
in understanding the magnetic properties of the (planar) 
doped Mott insulator in the stripe phase, and may serve 



as a microscopic basis to various spin models studied re- 
cently H[l0| • Although the model used by Zachar || dif- 
fers in many ways from the picture that we will present 
here, our conclusions arc in full agreement with his. 

It has been suggested both theoretically @Jll| and ex- 
perimentally Jl],||,^| that the charge ordering into hole 
rich stripes takes place before anti-phase spin domain 
modulations set in. In this spirit we will consider a 
model of the stripe phase of a doped AF Mott insulator 
in two dimensional square lattice based on the following 
assumptions: 1) the doped holes form an array of metal- 
lic stripes (domain walls) lining up, say in x-direction; 

2) between the charge stripes are AF insulating regions, 
each described by a Heisenberg spin ladder with rz leg legs 
of length N: 

N "log 

R S = J E E r) ■ [S(j + 1, r) + S(j, r + 1)] . (1) 
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3) the metallic stripes are Luttinger liquids in the Luther- 
Emery regime which has a spin gap. Thus, low energy 
hopping process between the stripe and the ladders is 
suppressed. || Our model system qualitatively corre- 
sponds to the electronic smectic phase of a doped Mott 
insulator first suggested by Kivelson, Fradkin and Emery 
|l2f , but we shall not include the effects of transverse fluc- 
tuations in the present paper. 

Now we apply the standard approach to a nonlinear 
cr-model (see, e.g., Ref. [Q) for the spin ladder. The 
partition function can be expressed as a path integral 
by means of the spin coherent states. For antiferromag- 
netic interaction (J > 0), the ladder is expected to be of 
short-ranged Neel order in the ground state. Therefore, 
we split the (staggered) spin fields into a slowly varying 
piece M, the order parameter field, and a small rapidly 
varying part I: 



S(j,r)~s (-y +r M(j) + al(j )r ) 
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where M ■ I = and M 2 — 1 . s is the spin (equal to 1/2 
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for the physical case) . We have assumed that the AF cor- 
relation length along the legs is much greater than the 
width of the ladder, which allows us to treat M constant 
along the rung s @. The low energy theory of the spin 
ladder is then found to be the familiar nonlinear er-modcl 
with the (Euclidean) action fig] 



dxdr 
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where the coupling constant and velocity are g a = 

\ («logErr' L rv) ^ C ° = aJs ( n log/ Err' L rr> ) ^ > 



and the matrix L r 
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The quantity 



"■("leg) = 0, 1 for ni og even and odd respectively. Ob- 
viously, the topological term appearing as the last in (||) 
vanishes for even-leg ladders which are in a Haldane gap 
phase 15 1. 

The charge stripe can be thought of a quasi-one dimen- 
sional electron gas (1DEG). In the intermediate temper- 
ature regime where the system is spin gapped but not 
superconducting yet, it becomes the Luther- Emery liq- 
uid with a (Euclidean) Lagrangian in the bosonized 
form 



-S^ideg — -§? c + -£s 



(4) 



The charge and spin pieces are each of the sine- Gordon 
variety: 
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where a = c, s for the charge and spin fields, respectively, 
and V c = (thus for cases of no Umklapp scattering). We 
consider only the Luther-Emery regime |l6f| of intermedi- 
ate temperature where Vs is relevant and a spin gap, A s , 
is dynamically generated. In this regime, the spin sine- 
Gordon theory is known to scale to a strong fixed point, 
K s = 1/2, sometimes called the Luther-Emery point. In 
the following, we shall restrict ourselves in performing 
calculations to this point, where some exact solution is 
possible. Our analysis can be perturbatively extended to 
around this fixed point. 

Now we are in a position to introduce the coupling 
between the charge stripe and the spin ladders on both 
sides (A,B): 
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^2[S A (j,ni eg ) + S B (j, 1)] ■ 40'>W(v0') (6) 



where f are the three Pauli matrices. Electrons in the 
stripe are coupled to the neighboring spins on the ladders 
at both sides. Such a spin-exchange interaction is quite 
natural from the study of the two dimensional t-J model 



for the high T c cuprates. Seeking for a low energy theory, 



we write c a (j) = \pla \e 
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where ipL.R are slowly varying left/right moving fermion 
fields with respect to the Fermi points ^kp, and bosonize 
them |l7| ]. The value of kp measures the filling factor of 
the stripe. We split each spin on the ladders into the 
slowly varying order parameter field M and the rapidly 
varying field I as in H mt splits into two separate 

parts, involving either M or I. We shall consider the 
part of M first and come back to another afterwards. 

The part of interaction involving the order parameter 
fields amounts to an extra term to the Lagrangian: 

■S?int,M = 9js(-) nies $ M (x) ■ 

where 4>m = Ma + VsMb for short-hand notation, and 
^2k F {x) are the spin piece of the usual 2kp spin den- 
sity wave order parameter (see, e.g., [H): 
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the dual field to 

vector. The relative phase factor 



^2k F = £MV2n9 s ), and = 
2w(j) s ) with a the short distance cutoff and 8 S 
Gjv = tt/cl is the Neel ordering wave 
77s, is +1 (—1) for site- 



centered (bond-centered) stripes, arising from staggering 
the spins of two neighbor ladders. The overall factor 
(_)™ic g appears simply because we have arbitrarily cho- 
sen the spin Sa (1, 1) as the unique reference vector when 
staggering spins for both ladders. As we shall see, this 
factor plays no physical role in the following calculation. 

Since the spin excitations in the stripe are gapped, the 
spin degree of freedom can be integrated out to obtain a 
correction term to the effective action of M; 
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To the lowest order in gj, 
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Here, the indices a,b = x,y,z. At the Luther-Emery 
point K s = 1/2, the sine-Gordon theory (||) for the spin 
sector is equivalent to a theory of free massive Dirac 
fermions [jl9|. The correlation function in (^|) can then 
be calculated exactly. At zero temperature, the result is 



(^(x,T)^ feF (0)) s = 



(10) 
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where K u (z) are Bessel functions. Here £ s is the spin 
correlation length related to the spin gap of the stripe, 
£ s = v s /A s , as introduced in Carlson et. al. ( |To| ) is 
manifestly spin rotationally invariant. 

Since the correlation function falls off at long distance 
(> £s) exponentially fast (due to the spin gap), we can 
use the operator product expansion to find an effective 
local interaction. The product of the (normal ordered) 
vertex operators e ±l "^ < t>" in (JsJ) at short distances yields 
a contribution proportional to the identity operator of 
the form 



g±iV27T0e (x,t) ^fiVZir&c (z 
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where we have dropped derivative operators. 

Inserting expressions ([To]) and ( JTT| ) to (||), one finds S' eS 
quadratic in 4>m (and so M's), in which the interaction 
falls off exponentially fast. Hence, the long wave length 
limit of (0) becomes (apart from an additive constant) 



S' a = 2Vrj s J dxdrM A ■ Mb 



(derivatives) (12) 



The effective coupling constant V is given by 
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The effective interaction term S' a is a correction to the 
nonlinear a- model described by the action S§ a in (^) . We 
should comment that we have used an effective long wave- 
length theory to determine the effective coupling Jl2] ) 
and that our expressions therefore depend strongly on 
the short distance cutoff. Also, a small parameter of 
Gn — 2kp cannot make our calculation more controllable 
for the same reason. This is natural since we are com- 
puting a non-universal quantity. 

So far we have left out the coupling between the 1DEG 
and the rapidly fluctuating fields I, which is part of the 
-ffint in (§)■ Analogous to ([?]), which comes from the in- 
teraction part of M, we find an interaction term j£?i n t,z 
of the fields I coupled to both the spin density and the 
2kp spin density waves. Nevertheless, the coupling to 
the latter becomes negligible for the rapidly oscillations 
in comparison with the former, so we ignore it. We then 
follow the same route as above to integrate out the spin 
degree of freedom in 1DEG, <f> s , and end up with effective 
couplings of I a and lg. Adding these to the terms ob- 
tained of each spin ladder itself, we find the total effective 
Lagrangian involving the fields I as follows 
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where the coupling constant V' — (gjs) 2 a/(4:TTV s ) is cal- 
culated exactly at the point K s = 1/2 from integrating 
out the spin fields 4> s of 1DEG. Since we generally expect 
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< 1 



we have reduced the problem of an infinite number of 
weekly coupled parallel spin ladders to a problem of two 



coupled ladders in (14) 
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FIG. 1. The (static) effective interaction of r\ a MA ■ Mb 
given by the integral expression (|13|). The coupling constant 



V is scaled with Vo = 




The lattice 
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spacing a is taken as the short distance cutoff such that inte- 
gration in the area of \/x 2 + v'ir 2 < a is excluded. Parameters 
are v c /v a = 3, and (a) K c = 1/4; (6) K c = 3/4. 

The effective Lagrangian is quadratic in the fluctuat- 
ing fields I, so they can be exactly integrated out, lead- 
ing to an effective coupling between the order parameter 
fields of two ladders, which is of the form 



X' ~ C"(M A x d T M A ) ■ (M B x d T M B ) 



(15) 



with C" a constant coefficient. Comparing with (|l2|), the 
above effective interaction is a higher order correction to 
the low energy effective theory. 

The static interaction of two adjacent spin ladders is 
solely given by (|l^) . The integral in ( |l3| ) has a short dis- 
tance divergence, which will be dealt with by imposing 
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a short distance cutoff. A natural cutoff is the lattice 
spacing a. Thus, our results will depend on the choice 
of cutoff, which is natural since this effective coupling is 
controlled by local physics. Obviously, additional oper- 
ator will yield corrections to the effective coupling but 
their effects should be considerably weaker. 

Fig. |l| shows that the effective interaction V changes 
sign when varying kp, namely the filling factor of the 
stripe, which for a free 1DEG are related by: kp = 
(1 — ph)ir/(2a) with ph the (linear) density of holes in 
each charge stripe. 

Table | summarizes what spin configurations are fa- 
vored by this effective interaction for different regions 
of kp . Our results suggest that the effective interaction 
could favor either in-phase or anti-phase spin modulation 
depending on whether the stripe is site-centered or bond- 
centered and on its physical parameters (such as line den- 
sity) . We would like to compare our results with others in 
two special cases. Firstly, we observe from Table | that 
for the bond-centered stripe the in-phase spin modula- 
tion is favored for small kp, i.e., high ph- Recently, Han, 
Wang and Lee studied the mean-field phases of the 
t- J+Coulomb model numerically. They found that for a 
doping range 0.02 < x < 0.14 the ground states are the 
bond- centered in-phase stripes, in which the line density 
of holes in each stripe varies from 0.38 at x = 0.095 to 
0.56 at x = 0.14, corresponding to a kp from 0.31^ to 
0.22-^ if assuming a free 1DEG. Since such a range of kp 
can be well considered below some kp, our result is con- 
sistent with that of Han, Wang and Lee (20) despite that 
their model and approach are completely different from 
ours. Secondly, we find for a site-centered stripe con- 
figuration that a transition from anti-phase to in-phase 
should occur at a critical kp (and thus a critical line 
density). That agrees with what Zachar |8| found in a 
different domain wall model, in which he estimated the 
critical (electron) filling factor to be 0.28 < 6 C < 0.30. 

The consistency of our result can also be checked by 
thinking of a special (but unphysical) case: half-filled 
(electron) charge stripe with kp = In this case, the 
2kp spin density wave of the stripe is commensurate with 
the short-ranged Neel order (with Gn = — ) of the neigh- 
boring AF spin ladders. Therefore, one expects that ef- 
fective coupling mediated by the stripe should become 
largest in this limit. Fig. [l] shows that the coupling con- 
stant V is indeed largest in magnitude for kp = ^L. 

TABLE I. Possible spin stripe modulations favored by the 
effective interaction (|l2]). Note that kp denotes qualitatively 
some phase transition point that can be read off from Fig. |l| 
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